Journal of Genetics and Genomics

Deep-sea corals are vital in maintaining coral ecosystem biodiversity, yet their genetic characteristics remain largely unexplored. Here, we present 11 deep-sea coral genome assemblies, including four Hexacorallia and seven Octocorallia species, significantly contributing new genomic information across two orders. Our analysis reveals the historical dynamics of coral speciation and the influence of environmental factors on the evolution of coral reef ecosystems.Total of 126 horizontal gene transfer (HGT) events were detected, among which genes from the ancestor of symbiodiniaceae indicate that the ancestors of deep-sea corals may have inhabited shallow-sea environments. Notably, several of these HGTs are involved in phosphorus (PhnX/PhnW) and cholesterol (DHCR7) metabolisms within corals, indicating that HGTs may serve as an adaptive survival strategy for the coral holobionts. Deep-sea corals also rely on symbiotic bacteria to synthesize 10 essential amino acids (such as valine and tyrosine), retaining only partial amino acid synthesis capacity. In addition, we investigated the evolution of key biological rhythm genes and temperature adaptation in corals. The loss of key rhythm genes (e.g., clock and cry) in deep-sea corals and copy number difference of genes related to heat stress (e.g., Cbl-b and Rchy) revealed genetic difference between deep-sea and shallow-sea corals. Our new genome assemblies enhance the understanding of deep-sea coral evolution, biodiversity, and adaptation, providing a genetic foundation for coral conservation.

WD40 domains share a widespread {beta}-propeller fold, and often act as versatile scaffold proteins. Despite their central role in organizing dynamic cellular complexes, the molecular and structural mechanisms of many WD40 proteins remain poorly understood. Among them, DCAF7, an ubiquitously expressed and essential gene in human, also encodes a highly conserved WD40 protein in eukaryotic organisms. It is known to interact with multiple and functionnally diverse partners to coordinates cellular activity of several protein kinases as well as transcriptional regulators, thereby modulating key cellular processes such as cell growth, differentiation, and transcriptional regulation. However, the precise mode of action of DCAF7 is unknown and its important divergence in sequence from better characterize WD40 prevent information transfer by similarity. Structural interactomic can reveal how protein-protein interactions (PPIs) occur within an organism and are essential for understanding biological functions and developing new therapeutic strategies. Using SLiMAn2, AlphaFold2/3 and PSSMsearch, we identified a conserved -helical short linear motif (SLiM) in several well known DCAF7 partners that binds to the top surface of its {beta}-propeller. This motif was subsequently used to generate a regular expression, to identify potential new direct binders across the DCAF7 meta-interactome and the human proteome. Domain-domain interactions were also predicted for some other partners. Finally, modeling of oligomeric complexes with such new hits reveals the structural basis of DCAF7 scaffolding, with links to neurodevelopmental disorders such as autism.

Achieving high-throughput and precise phenotypic quantification and imaging modalities of stomatal and epidermal cells across diverse species remains a primary bottleneck in elucidating the mechanisms of stomatal dynamics, epidermal patterning, and environmental adaptation of plants. Here, we developed EpiReasoner, an artificial intelligence framework comprising a vision module, EpiVision, and a knowledge-based reasoning module, EpiBrain, for the quantitative phenotypic analysis and domain-specific knowledge reasoning of stomatal complexes and pavement cells in plants. Operating across bright-field, scanning electron microscopy, and differential interference contrast modalities, EpiVision achieves precise instance segmentation in various monocotyledonous, dicotyledonous, and fern species. Its performance significantly surpasses current state-of-the-art models. Moreover, we defined 23 quantitative indices describing stomatal cell morphology and spatial distribution. For domain-specific tasks such as phenotype prediction, genotype deduction, and molecular mechanism reasoning, EpiBrain demonstrates a human preference rate significantly higher than that of general-purpose large language models, including GPT-5 and Claude Sonnet 4. The application of EpiReasoner to phenotypic data of stomatal density derived from a tomato natural population of 170 accessions successfully identified a major quantitative trait locus on chromosome 8. The candidate gene, SKP1-interaction partner 19L (SKIP19L), encoding an F-box family protein, exhibited severe allele frequency drift during tomato domestication, which is highly consistent with the adaptive trend of reduced stomatal density under artificial selection. EpiReasoner provides a novel paradigm that unifies visual phenomics and knowledge-driven reasoning for the biology of stomata and pavement cells, thereby significantly accelerating scientific discovery in plant science.

The Caenorhabditis elegans AWC olfactory neuron pair specifies asymmetric subtypes, AWCOFF and AWCON, through stochastic and coordinated cell signaling events. UNC-104/kinesin-3 (KIF1A) and UNC-116/kinesin-1 motor proteins act in the AWCON cell to regulate the synaptic localization of the TIR-1/SARM1-assembled calcium signaling complex in the AWCOFF cell to promote AWCOFF. However, the molecular mechanism in the AWCON cell that acts non-cell autonomously to control synaptic TIR-1 calcium signaling to promote AWCOFF remains unclear. Here, we show that JIP-1, a conserved c-Jun N-terminal kinase (JNK)-interacting protein 1, mediates the synaptic localization of TIR-1 in the AWC axon to specify the AWCOFF subtype. A jip-1 loss-of-function mutant, identified from an unbiased forward genetic screen, has reduced localization of TIR-1 at synapses in the AWC axon and accumulation of TIR-1 in the AWC cell body. jip-1 mutants significantly enhance the 2AWCON phenotype of a hypomorphic tir-1 mutant. JIP-1, like UNC-104 and UNC-116, mainly acts non-cell autonomously in AWCON to specify the AWCOFF subtype. Our findings provide mechanistic insights into how cell-specific Ca2+ signaling proteins, such as TIR-1, target synaptic regions via intercellular signaling to promote neuronal diversification.

Centrosome dysfunction is linked to developmental disorders affecting brain and body size, including microcephaly and primordial dwarfism. However, the cellular mechanisms underlying these rare conditions remain poorly understood. In this study, we investigate a rare variant of the centrosome-associated protein Pericentrin, which was discovered in a single family with Majewski/microcephalic osteodysplastic primordial dwarfism type II (MOPD II). Unlike the majority of pathogenic PCNT variants that cause severe protein truncation, the p.Lys3154del variant ({Delta}K3154) involves a single amino acid deletion in the proteins only conserved functional domain, providing a unique opportunity to explore PCNT function in MOPD II. To model PCNT{Delta}K3154, we examined the effects of Drosophila Pericentrin-like protein (PLP) carrying an orthologous deletion (Plp{Delta}R). Our results show that plp{Delta}R animals exhibit smaller tissues that recapitulate MOPD II phenotypes. Behavioral assays revealed defects in climbing and mechanosensation, suggesting impaired sensory cilia function. We also found that Plp{Delta}R cells exhibit accelerated mitosis, increased apoptosis, and reduced pericentriolar material recruitment. In silico structural modeling, yeast two-hybrid, and co-immunoprecipitation experiments show that Plp{Delta}R produces a protein that disrupts PLP dimerization and PLP interaction with Asterless, another centrosome protein. Overall, modeling the human MOPD II patient variant PCNT{Delta}K3154 in Drosophila reveals how a single amino acid deletion affects biological processes from the molecular level to the organismal level. Our work offers new insights into the defective cellular mechanisms underlying MOPD II in patients with the PCNT{Delta}K3154 variant, potentially linking the etiology of the disease in these individuals to the loss of a single protein-protein interaction.

Heterogeneous nuclear ribonucleoprotein U (HNRNPU) deficiency is a rare genetic cause of neurodevelopmental disorders (NDDs) lacking targeted therapies. Here, we developed a transcriptomic-guided compound prioritization pipeline using Connectivity Map (CMap) analysis on multi-model transcriptomic signatures from HNRNPU-deficient human cells and mouse models. Ten compounds were selected through manual curation and functionally screened in patient-derived HNRNPU-deficient neuroepithelial stem (NES) cells with earlier observed cellular phenotypes. Two of the compounds, AS601245 and Lenalidomide, significantly reduced the elevated neural progenitor population during differentiation, and their combination further decreased primary cilia incidence, indicating partial rescue of the patient-specific cellular phenotypes. To understand the mechanisms underlying the partial rescue, we employed proteome integral solubility alteration (PISA) and expression proteomics. PISA assay identified TMEM150C and GSK3A as proximal targets of combined treatment. Additionally, we observed reversal of multiple biological pathways including downregulation of Wnt signalling and upregulation of mitochondrial pathways and transmembrane proteins. Altogether, we established a computational-experimental pipeline for transcriptomic-guided drug repurposing for a monogenic NDD, and demonstrated that the network-level modulation partially rescues the delayed neural differentiation in HNRNPU-deficient neural cells.

Climate change poses an increasing threat to the cultivation of deciduous fruit trees, placing greater demands on modern pear breeding. Using pear germplasm adapted to diverse environments, we assembled 11 chromosome-level genomes. In combination with 13 publicly accessible pear genomes, we analyzed presence-absence variations (PAVs) and constructed a graph-based pangenome for pear. By performing a PAV-eQTL analysis of the fruit of 123 pear accessions, we identified PAVs significantly associated with expression levels of genes that may be involved in regulating agronomic traits. Population analysis of 268 pear accessions revealed two stop-gained variants in DAM1 of independent origin, which may function in advancing the blooming date and reducing the chilling requirement. We detected complex PAVs at the NOR1 locus, including two copy-number variations and one deletion. These PAVs contributed to the rapid diversification of the NOR1 locus and the fruit development period through regulating ARF5 and other ripening-related genes. We revealed the selection history of the NOR1 locus and developed novel pear individuals that accumulated alleles for low chilling requirement, early blooming date, and short fruit development period. The results provide valuable resources for pear genomics research and offer a guideline for breeding modern pears with climate resilience.

Abstract: Objective This study aimed to systematically screen for potential candidate biomarkers and identify therapeutic targets associated with gouty arthritis (GA) through integrated analyses of single-cell and bulk RNA sequencing (RNA-seq) data. Methods The single-cell dataset GSE211783 and the bulk RNA-seq dataset GSE160170 were analyzed using a series of bioinformatic approaches, including cell clustering, differential expression analysis, immune cell infiltration assessment, protein-protein interaction network construction, gene set enrichment analysis, as well as drug sensitivity evaluation. To establish an animal model of GA, monosodium urate crystals were injected intra-articularly into experimental mice. Joint swelling was evaluated, and morphological changes in joint tissues were analyzed through hematoxylin-eosin staining. The presence of TREM1-positive cells was detected by immunohistochemistry and the level of TREM1 protein expression in joint tissues were assessed by Western blotting. Results We identified 102 differentially expressed genes (DEGs) and 14 signaling pathways associated with GA. The PPI network revealed 25 hub genes, of which 17 (including TREM1, TNF, PTGS2, and NLRP3) were highly expressed and 8 (including FCGR3B and CXCR6) showed low expression in the GA samples. These genes correlated significantly with the infiltration levels of macrophages. Among the hub genes, TREM1 was selected for further validation because it correlated significantly with all 14 differential pathways. In animal experiments, GA mice developed marked joint swelling and inflammatory tissue injury, along with a significant increase in TREM1-positive cells and TREM1 protein expression. Conclusion Integrative analysis of single-cell and bulk RNA-seq data identified 102 GA-related DEGs and 14 key pathways, from which 25 hub genes were screened. TREM1 is significantly upregulated in GA and may be linked to macrophage function, providing new insights into biomarker and therapeutic target discovery for GA.

The Japanese quail (Coturnix japonica) is a widely used model organism in developmental biology, genetics, and agriculture. Here, we present new, haplotyped, high-quality genome assemblies of the Japanese quail, generated using a combination of state-of-the-art sequencing technologies, including PacBio HiFi long reads, Oxford Nanopore sequencing, and Hi-C scaffolding. This assembly has a total length of 1.19 Gb, 80% of which is included in chromosomes, and is highly complete (BUSCO score aves_odb10: 97.3). Assembly metrics show a marked improvement in contiguity, with a significantly higher scaffold N50 and a lower number of contigs compared to the reference genome assembly. Remarkably, the assembly extends previously truncated chromosome ends, with 31 telomeres detected. In addition, we merged the existing Ensembl and Refseq annotations and obtained a combined set of 26,102 genes, of which 25,038 genes were successfully mapped on the improved assembly haplotype 1 (Cjap1.hap1). Together, these new genome assemblies and their enriched annotation provide a robust genomic framework for future research. They enhance our ability to investigate developmental processes, genetic and epigenetic inheritance, and host-pathogen interactions. Furthermore, they offer valuable insights for conservation genetics and sustainable breeding programs. This resource represents a critical step forward in leveraging the full potential of the Japanese quail as a model species in both basic and applied research.

PUF60 is a splicing factor related to the polypyrimidine-tract binding protein U2AF2. PUF60 is deleted in developmental disorders such as Verheij syndrome and amplified in approximately 8% of cancers. Thus, both increases and decreases in PUF60 expression can have profound physiological effects. However, little is known about how changes in PUF60 expression impact global splicing patterns. Here, we created a model system of CRISPRa/i in mouse stem cells (mESCs) to transcriptionally upregulate or downregulate Puf60. Our results uncovered extensive transcriptional, post-transcriptional, and post-translational regulation of Puf60 protein expression. We observed that Puf60 protein levels in normal mESCs drop dramatically at a critical cell density, leading to cell death. Puf60 is very essential in stem cells, and its repression causes cell death and impacts specific splicing events, including its own splicing autoregulation, providing valuable insights into the functional consequences of PUF60 dysregulation. Analysis of phosphoprotein data revealed phosphorylation of threonine at the N-terminus of PUF60. Our results showed that mutating threonine to glutamate downregulates the protein and alters its localization. Thus, our study reveals a novel regulatory mechanism of Puf60 phosphorylation that mediates its function and may be related to its frequent overexpression in cancer cells.

Human papillomavirus 18 (HPV18) preferentially infects cervical stem cell-like cells and is strongly associated with adenocarcinoma. However, the mechanisms underlying differentiation into cervical adenocarcinoma remain unclear due to the lack of appropriate experimental models. We aimed to establish a model of HPV18-associated cervical adenocarcinoma and elucidate its molecular and cellular differentiation mechanisms. HPV18 E6/E7 were introduced into induced pluripotent stem cell-derived reserve cell-like cells (iRCs) to generate tumor models. Spatial transcriptomics and single-cell multi-omics analyses were performed to integrate histological and molecular data. A distinct component (Gland_A) exhibited morphological and immunohistochemical features of cervical adenocarcinoma and was efficiently induced in iRC-18 tumors. Gland_A showed increased chromatin accessibility and elevated expression of FOXA1, FOXA2, and ALDH1A1. Analysis of clinical samples confirmed enrichment of ALDH1A1 in HPV-associated adenocarcinomas. This model recapitulates key features of HPV18-associated cervical adenocarcinoma and provides insights into its differentiation mechanisms.

Magnaporthe oryzae, the rice blast fungus, plays a role as a model organism for molecular plant-microbe interaction research. Studies on the pathogenic mechanism of this fungus revealed many genes involved in signaling pathways. As multi-omics data are being available, genomic-level researches have been conducted to uncover the underlying biological processes during the pathogenesis of M. oryzae. Identifying the genome-wide protein-protein interaction (PPI) network is one of the omics-level approaches, which helps to understand signaling and regulatory pathways. However, existing biological network resources of M. oryzae are not sufficient to decipher pathogenesis mechanisms due to the abundance of false positives/negatives. In this study, a reliable PPI network database of M. oryzae, MagNet, was constructed with three methods, including homology-based Interolog search, co-expression network construction, and domain-domain interaction (DDI)-based prediction. With three approaches altogether, the pan-network with 5,600,976 interactions was generated, including 217,531 highly confident interactions supported by all three methods. Experimental data on M. oryzae PPIs supported that our PPI network can predict PPIs with higher accuracy compared to the previously constructed databases. MagNet would provide integrated biological network data, which can help to understand the molecular mechanisms of the rice blast fungus. The PPI data can be accessed via https:/magnet.scnu.ac.kr.

In wheat, weak seed dormancy (SD) is related to an increased tendency for pre-harvest sprouting (PHS), which reduces yield and quality. However, the molecular mechanism underlying SD remains elusive. Here, we identified a wheat R2R3-MYB transcription factor (TaMYB83-7B) related to SD. Expression analysis showed that TaMYB83-7B was highly expressed in wheat seeds, and was more highly expressed in strong-dormancy varieties than in weak-dormancy varieties. Sequence and association analysis indicated that T/C mutations at -907 bp and -1133 bp in the TaMYB83-7B promoter were significantly associated with wheat SD, with C at both sites related to strong dormancy. Dual-luciferase reporter assays demonstrated that the transcriptional activity of the TaMYB83-7B promoter was significantly higher in strong-dormancy varieties than in weak-dormancy varieties. Further analyses indicated that TaMYB83-7B functions as a transcriptional inhibitor. Germination experiments revealed that overexpression of TaMYB83-7B significantly enhanced SD, while its loss-of-function reduced SD. Finally, TaMYB83-7B was found to regulate SD by influencing the balance between abscisic acid (ABA) and gibberellin (GA) in wheat seeds. Overall, the results of this study enhance our understanding of the complex regulatory mechanism underlying SD, and provide gene targets and molecular markers for the genetic improvement of PHS resistance in wheat.

Plants continuously develop shoot branches derived from axillary meristems. In rice (Oryza sativa), TILLERS ABSENT1 (TAB1), an ortholog of Arabidopsis WUSCHEL, plays an essential role in axillary meristem formation by promoting stem cell proliferation. Although several genes associated with TAB1 function have been identified, the molecular mechanisms underlying stem cell proliferation during axillary meristem formation remain poorly understood. Here we identify ABERRANT SPIKELET AND PANICLE1 (ASP1), a TOPLESS-like transcriptional corepressor, as a novel regulator of axillary meristem formation, and investigate downstream mechanisms regulated by TAB1 and ASP1. In asp1, the stem cell region was expanded, indicating that ASP1 negatively regulates stem cell proliferation. Notably, WOX4, a paralog of TAB1, was precociously expressed in asp1, possibly in association with expansion of the stem cell region. Genetic analysis further revealed that asp1 mutation rescued the loss of axillary meristems in tab1. Transcriptome analysis showed that several type-A RESPONSE REGULATOR (OsRR) genes, encoding negative regulators of cytokinin signaling, were upregulated in tab1 relative to wild type, asp1, and the tab1 asp1 double mutant. Consistently, fluorescence of the synthetic cytokinin reporter was absent during axillary meristem formation in tab1 but was detected in wild type and tab1 asp1. Moreover, overexpression of OsRR10 inhibited axillary meristem formation, phenocopying tab1. Collectively, these findings suggest that TAB1 activates cytokinin signaling by repressing type-A OsRR expression, whereas ASP1 negatively regulates cytokinin signaling by promoting the expression of these genes. Thus, rescue of the tab1 phenotype by asp1 mutation probably reflects restoration of cytokinin signaling.

O_LISymbiosis between legumes and rhizobia is beneficial on nutrient-poor soils, as it enables the fixation of atmospheric N2. To establish this symbiosis, gene expression in both the host plant and the symbiont has to be regulated. To understand the underlying RNA-mediated regulation of host gene expression, we designed experiments to identify competing endogenous networks involving circular RNA, microRNA, and linear transcripts during symbiosis, using wt and symbiosis-deficient Lotus japonicus mutants with the rhizobium Mesorhizobium loti (M. loti). C_LIO_LICircRNA, miRNA, and linear transcripts were identified from Lotus japonicus wildtype and CCamK mutant (ccamk-13; snf-1) seedlings without inoculation or with M. loti inoculation using deep short-read sequencing with rRNA-depletion and random primers. C_LIO_LIDifferentially expressed miRNAs showed negative correlations to predicted target genes and may regulate symbiotic processes. The symbiosis essential iron-sensor LjnsRING/BRUTUS expresses a circRNA which was upregulated in symbiotic treatments. This circRNA may act as a target mimic and contribute to nodule longevity. CircRNAs are predicted to act predominantly as trans-regulatory molecules with similar frequencies in Arabidopsis thaliania, Oryza sativa, and Lotus japonicus. C_LIO_LIWe identified novel miRNAs, long noncoding RNAs, and circRNAs, and nominated several as potential new regulatory non-coding RNAs that may act as target mimics to stabilize genes and support symbiosis. C_LI SummarySymbiosis between Lotus japonicus and Mesorhizobium loti involves treatment-specific regulation of competing endogenous RNA networks involving circular RNA, miRNA, and linear transcripts.

Histone lactylation is a recently identified histone post-translational modification (PTM) that links energy metabolism to chromatin regulation. Although histone lactylation has been implicated in transcriptional activation, its function in meiotic chromatin remains unclear. Previously, we identified enrichment of multiple histone lactylation marks within the meiotic karyosome, a highly condensed and transcriptionally repressive chromatin structure formed in Drosophila oocytes. Here, through an RNAi-based screen, we identified the CBP family protein dCBP as a regulator of histone lactylation in the karyosome. Germline-specific knockdown of dCBP preferentially reduced histone lactylation, including H4K8 lactylation, and caused premature disruption of the synaptonemal complex, abnormal egg chamber development with excess nurse cells, reduced egg production, and decreased embryonic viability. Corresponding histone acetylation marks were comparatively less affected than histone lactylation by dCBP knockdown. Together, our findings provide evidence that dCBP-mediated histone lactylation contributes to meiotic chromosome maintenance and suggest a potential link between energy metabolism and meiotic chromatin regulation.

Background: Inherited Retinal Diseases (IRDs) are a group of genetically heterogeneous blinding conditions. Major global genomic reference databases are disproportionately enriched for individuals of European ancestry. This underrepresentation creates a significant bias that impedes the accuracy of genetic diagnosis in the Chinese population. This study aims to address this limitation by constructing a comprehensive genetic landscape of IRDs using large-scale deep-sequencing data from a large Chinese cohort. Methods: The study leveraged variant data primarily from 10,588 individuals in the China Metabolic Analytics Project (ChinaMAP) and cross-referenced findings against multiple national and international databases. We systematically curated variants within a targeted panel of 291 IRD-associated genes. Variant pathogenicity was assessed using a comprehensive pipeline integrating InterVar-automated classification based on 2015 American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines, ClinVar evidence (review status [≥] 1 star), and manual literature curation. We delineated the mutational spectrum, identified population-enriched pathogenic/likely pathogenic (P/LP) variants, and analyzed the distribution characteristics of IRD-associated highly-mutated genes. Furthermore, we calculated the carrier frequencies (CF) and genetic prevalence (GP) of autosomal recessive(AR)-IRD genes in the Chinese population. Results: The study revealed a highly concentrated genetic landscape for AR-IRDs in the Chinese population, with ABCA4 and USH2A emerging as the primary drivers of the genetic burden. This finding aligns with previous Chinese cohorts but contrasts with global databases, where genes such as the X-linked RPGR are more prevalent. In contrast, autosomal dominant (AD)-IRDs exhibited high locus heterogeneity, with pathogenic variants dispersed across numerous genes (e.g., COL2A1 and MFN2). We identified a series of P/LP variants that were either high-frequency or significantly enriched in the Chinese population, such as CNGB1 (p.P530R) and specific recurrent alleles in ABCA4 and CYP4V2. The estimated cumulative CF for AR-IRDs was 1 in 5.60, and the theoretical total GP was 1 in 2,624.67, based on the ChinaMAP data. Conclusion: By integrating the ChinaMAP dataset with diverse genomic resources, this study provides a genetic landscape of IRDs in the Chinese population. Our analysis shows a concentrated mutational spectrum in AR-IRDs, contrasting with the pronounced heterogeneity in AD-IRDs. These findings, including population-specific pathogenic variants and refined prevalence estimates, provide a resource for precision diagnostics, genetic counseling, expanded carrier screening (ECS), and public health policy development in China.

Accurate modeling of genotype-by-environment (GxE) interaction is critical for genomic prediction in plant breeding but remains challenging due to complex interaction structures. Conventional models often use the Hadamard product of genotype and environment covariance matrices to capture joint similarity, which may not fully represent GxE complexity. Here we propose a novel framework that derives covariance structures from the matrix multiplication of genotype and environment kernels, decomposing these into symmetric components incorporated as random effects in mixed models. Evaluated for 11 wheat and rice multi-environment datasets and across, this approach consistently outperformed the traditional Hadamard-based model, improving prediction accuracy by up to 13.2% in Pearsons correlation and enhancing top-selection accuracy. Combining both methods yielded the highest performance, indicating complementary information capture. This framework offers a flexible, interpretable, and computationally feasible extension for modeling GxE interaction, potentially enhancing genomic selection effectiveness under diverse environmental conditions.

Copy number variations (CNVs) are major structural genomic variants that contribute to a wide range of human diseases. Accurate detection of CNVs from whole-exome sequencing (WES) data has been a long-sought goal for clinical and population genetic studies. Despite recent progress, existing WES-based CNV callers still suffer from high false-positive rates and reduced recall for short-length variants, and current deep learning methods have not fully used complementary information in region-level genomic features. Here we present CN-RNN, a deep learning-based CNV caller for WES data. The model combines a bidirectional long short-term memory (BiLSTM) branch that captures local depth changes and contextual dependencies across neighboring exons with a parallel multi-layer perceptron (MLP) branch that encodes region-level metadata such as GC content, mappability, and exon length. CN-RNN was trained on the Autism Sequencing Consortium (ASC) parent-child trio cohort using the Mendelian rule of inheritance to ensure high-quality training sets. It was evaluated across three independent datasets, in which we showed that CN-RNN outperformed existing WES-based CNV callers and deep learning methods. CN-RNN offers a scalable, accurate tool for CNV profiling in WES-based studies and supports broader application of CNV analysis in population and clinical research. CN-RNN is available at https://github.com/FeifeiXiao-lab/CN-RNN.

The ability to generate functional B cells from human pluripotent stem cells (hPSCs) would open new opportunities to develop novel B cell-based therapies to treat a range of human diseases and disorders. Towards this goal, we established a protocol that promotes the efficient development of B lineage cells from definitive hematopoietic progenitors generated from different hPSC lines. Flow cytometric and multi-omic scRNA-seq analyses revealed that B cell development from hPSCs transitions through the well-established pro-B, pre-B and naive B cell stages, accurately recapitulating B lymphopoiesis in the human adult bone marrow. Importantly, the naive B cells generated with this approach could be induced to mature into plasma cells that secrete antibodies and undergo class switching. Analyses of signaling pathways that regulate B lymphopoiesis in these cultures uncovered a potent inhibitory effect of IL-7 on functional IgH rearrangement, resulting in the development of abnormal cells that failed to undergo pre-B cell maturation. Finally, analysis of the different hPSC-derived hematopoietic programs revealed that both definitive and yolk sac progenitors display B cell potential, indicating that there are distinct developmental sources of human B lineage cells. Taken together, these findings demonstrate the efficient generation of B cells from hPSCs and, in doing so, provide a system for further investigating the earliest stages of human B lymphopoiesis and a source of appropriately staged plasma cells for future therapeutic applications.